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ABSTRACT 

Cosmic shear measurements rely on our ability to measure and correct the Point 
Spread Function (PSF) of the observations. This PSF is measured using stars in the 
field, which give a noisy measure at random points in the field. Using Wiener filter- 
ing, we show how errors in this PSF correction process propagate into shear power 
spectrum errors. This allows us to test future space-based missions, such as Euclid or 
JDEM, thereby allowing us to set clear engineering specifications on PSF variability. 
For ground-based surveys, where the variability of the PSF is dominated by the envi- 
ronment, we briefiy discuss how our approach can also be used to study the potential 
of mitigation techniques such as correlating galaxy shapes in different exposures. To 
illustrate our approach we show that for a Euclid-like survey to be statistics limited, 
an initial pre-correction PSF ellipticity power spectrum, with a power-law slope of -3 
must have an amplitude at £ — 1000 of less than 2 x 10~^^. This is 1500 times smaller 
than the typical lensing signal at this scale. We also find that the power spectrum of 
PSF size ((5ij2) at this scale must be below 2 x 10^^^. 
Public code available as part of iCosmo at http://www.icosmo.org 
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1 INTRODUCTION 

High accuracy cosmic shear measurements require a precise 
measure of galaxy shapes, which become correlated due to 
gravitational lensing. To measure this lensing induced corre- 
lation, we must first correct for any instrument effects that 
also cause the observed galaxy shapes to correlate. In par- 
ticular, the Point Spread function (PSF) of the instrument 
needs to be corrected. This PSF includes effects associated 
with the instrument as well as the observing environment 
such as the atmosphere. Given an observing system with a 
given PSF, correcting a galaxy image is done by measur- 
ing the PSF using stars in the neighborhood of the galaxy 
and then deconvolving this from the galaxy image. In so do- 
ing, the two point correlation function of galaxy ellipticities 
(which is a combination of true galaxy ellipticity and PSF el- 
lipticity correlation functions) can be reduced to true galaxy 
ellipti city with some residual error. In|A mara and Refregieij 
(|2008l ) we showed that future ambitious surveys will need 
to control these residuals to a level where their contribution 
variance over the range of scales being used is a^y^ < 10~^, 
where a'^y^ is defined in lAmara and R cfrcgicr (2008). This 



results was also confirmed by iKitching et al.l (|2008l ') 



Achieving this will require both tight controls of the in- 
herent PSF correlations of the instr ument as well as an accu- 
rate method for PSF correction. In lPaulin-Henriksson et abl 



(|2008l ) we focused exclusively on the limits of PSF correc- 
tion without considering the impact of the initial PSF pat- 
tern. This allowed us to investigate the minimal amount of 
information that we would need to collect in order to cor- 
rect the PSF so that the residuals stay below crj^^ < 10~^. 
Since each star gives a finite amount of information about 
the PSF due to noise and pixelisation, the PSF must be 
measured by collecting information from a number of stars. 
The length scale that corresponds to this minimal number 
of stars ( Omin) then becomes a key scale, where correlation 
information of larger scales is 'safe' but ellipticity correla- 
tions on scales smaller than this are 'not safe' unless ex- 
tra information is known a bout the PSF. In this pa p er, we 
extend the earlier work in iPaulin-Henriksson et al.l (|2008| ) 
, using Wiener filtering, to include the requirements that 
would need to be placed on the PSF correlation function 
if we wish to use scales smaller than Omin- In this way the 
design of future surveys can be divided into two regimes: (i) 
large scales where the adverse affects of PSF are correct in 
software using image processing; and (ii) small scales that 
need to be controlled in hardware by, for instance, designing 
a space-based instrument with a well behaved PSF on small 
scales. 

This paper is organised as follows: Section 2 sets out 
the basic statistical properties of a field that is sampled at 
random points. In section 3, we show the impact of using a 
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Wiener filter to model the PSF which leaves residual errors. 
In section 4, we show the impact for cosmic shear surveys, 
and finally we discuss our interpretation and conclusions in 
section 5. 



^}f = {\fm') = j^ j MC^n^), (8) 
This can be substituted into Equation [7] to give 



(9) 



2 RANDOM SAMPLING 

To correct an image for the effects of the PSF we must be 
able to construct a model for the PSF at the desired point, 
which is usually the position of a galaxy. We do this by 
measuring the PSF, which we assume is a continuos func- 
tion f{0), at discrete points Gi where the stars lie in the 
image. We can construct a representation of the observed 

field fobsiey. 



(1) 



where Tj{9) is a sampling operator that keeps information 
only at the position Gi, given by 



Ns 



■0^). 



(2) 



Here S is the Dirac delta function and Ns is the total number 
of stars over an area A. Equation [T] also contains measure- 
ment errors - e, with (e) — and (e^) = ct^^, which is also 
sampled at the positions of the stars. With a representation 
of the observed field, we can now study the statistical prop- 
erties of the field, namely the two-point correlation function 
of a field /, 



m = {f{e)f{e 



M \fi£)\' 



-i£.it> 



(3) 



which is the well known result that the correlation function 
is the Fourier transform of the power spectrum C^^{£) = 
!/(.«) p. From Equation [U we can define the functions g{0) = 
f{0)rj{9) and h{9) — e{9)rj{9). This allows us to decompose 
the observed power spectrum, C°'"'(^), into three compo- 
nents, C»»(^), C''''(£) and C<'''[t) that are, respectively, the 
power spectrum from the autocorrelation of g{6), the au- 
tocorrelation of h(6) and the crosscorrelation of g{6) and 

HO), 



(4) 



The power spectrum C^^ is constructed from the Fourier 
space expression of g{9), 



m 



(2^)2 
using the fact that 



d£'f(£')fj{£-£') 



V{£) 



A ^= 
i = l 



The power spectrum from autocorrelation of g{d) is 



C'"'(£) = C^^(£) + 



n,(27r)2 



d£'C^^{£'), 



(5) 



(6) 



(7) 



where Ua = Ns/ A is the star density. The integral in Equa- 
tion [7] can be expressed in terms of the point variance of 



Applying the same procedure to the rest of the terms in 
Equation!?] leads to the expression, 

2 2 2 

C°'"'{£) = C"{£)+Ci£)+2C'''{£) + ^ + ^+2^.{W) 

We see that our observed correlation function is com- 
posed of three terms that come from the correlations and 
cross-correlations of the underlying field of interest and the 
errors, and three terms that look like white noise (scale in- 
dependent) terms. 



3 WIENER FILTERING 

The expression shown in Equation [10] gives the raw power 
spectrum that would be measured. In practice, we can con- 
struct an estimator of the field / that is better behaved. 
Here we will do this using Wiener filtering, which is known 
to be an optimal filter for reducing the for the residuals. 
The estimator / can be construct as 



fi£)^mfoU£), 

where the filter $ is given by 



(11) 



(12) 



with the observable being defined as C"*" {£) = C^^ {£) + 
C"{£) and C"'{i) is the power spectrum of the noise (see 
equation [To} . From this we can calculate the residuals, 

slw = ?{£) - m = musW - m, (13) 

and the power spectrum of the residuals, 

C^f»f(^i) = <1>^(^)C°'"(£) + C"{£) - 2^£)C°'"^{£). (14) 

If we assume that the measurement error, e{£), is not 
correlated to the signal, f{£), in Equation[l]then the term for 
the observed power spectrum in Equation[TD]can be reduced 
to. 



C°'"{£) 



C^^{£)+C{£) + ^ + ^ 
ris Us 

C"{£)+C"{£), 



(15) 
(16) 



where C"(£) contains the noise terms, which in this case are 
C"{£) = C''{£) + crf f/ns + ol^lus. In this case the power 
spectrum of the residual in Equation [14] reduces to 

CW(^) = ($(^) - l)'c-^^'(^) + $2(^)C"(£). (17) 

For simplicity, in what follows we will consider the isotropic 
case where there is no preferred direction. In this case, the 
powerspectra become dependent only on i and so C{€) = 
C{£). In this isotropic case, we can substitute the definition 
of the Wiener filter in Equation [12] into Equation 1171 which 
leads to 
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(18) 



4 LINK WITH COSMIC SHEAR 

In IPaulin-Henriksson et al] (|2008l ) we showed how errors in 
the measurement of PSF size and ellipticity propagate into 
errors on the measured shear. This can be expressed as 

Sf^FSj,2+GS„ (19) 

where 3^2 is the variation of the PSF size, 5^2 — {Rp — 
(Rp)) / (Rp) and similarly Se is the variation in PSF elliptic- 
ity. F and G are given by: 



Rp 

V Rn 



(Cg Ep), 



and 



_ 1 f Rp 

^-~p^[r: 



(20) 



(21) 



Here, P~' is the shear susceptibility factor, R and e are the 
Radius and ellipticity and the subscripts p and g stand for 
PSF and galaxy, respectively. 

For the rest of this paper, we will focus on the case 
where on average the errors in the two component of shear 
(71 and 72 where 7 = 71 + 172) are the same, and thus 
the power spectrum can be written in terms of one of the 
components, 

*7 2C^^^ 



(22) 



One of the shear components (for instance the first one) is 
given by, 

5^1^ Fi5r2+G5,^, (23) 

which leads to a systematic power spectrum of 

Cp^ = Fld'^'' + G^Cp^ . (24) 

Both R^ and ei are fields that are constructed from the 
measurements of the stars. From Equation ll8l we know that 
the post Wiener filtering power spectrum of these two fields 
is given by 



c, 



(25) 



where y can be eithe r R? o r ei. The systematic limit from 
lAmara and Refreeierl ()2008l ) can then be expressed as 



2at 



(26) 



(27) 



4.1 Illustrative Examples 

We now explore some example cases to illustrate how an 
analysis of the PSF variation can be set by a desired sys- 
tematic target. For this we will make some simplifying as- 
sumptions. First, we will assume that the systematic contri- 
butions from PSF size are comparable to those coming from 
ellipticity. Explicitly, we set o-^ygji2 = <^1ysti ■ We also set 



P^ = 1.84, (eg)^/^ = 0.4 {eDy ^ = 0.05 and Rg/Rp = 1.5 
(iPaulin-Henriksson et al. I I2OO9I ). We also assume that the 
star density if 1 star per arcmin^ with a S/N = 100, 
which in turn can be translated to errors on the PSF, per 
star, of o-(e) = 3 x 10"^ and a{R^)/R^ = 3 x 10~^ (see 
IPaulin-Henriksson et al.l (|2009l l for justification) . These are 
conservative assumptions for a Euclid like mission, but we 
should should note the our requirements will depend, mod- 
erately, on the choice of these parameters. Given this in- 
formation, we can set requirements on the powerspectra of 
the PSF variation that are needed to reach a systematics 
level of alys < 10"'^, given a particular functional form. For 
simplicity we begin with a power-law 




(28) 



where to is some convenient reference scale, we pick lo — 
1000. Note that the choice of £0 has no impact on the results, 
for instance those shown in figure[T] This choice only chances 
the Gf that we quote. Gf is the amplitude of the power 
spectrum at this £ value and /3 is the slope of the power-law. 
In Figure [T] we see that the most relaxed requirements for 
Gf come at roughly P — —3. This is because for slopes shal- 
lower than this, the contribution to that variance from small 
scales (high £) begins to dominate, whereas for steep spec- 
tra the variations on large scales start to become a problem. 
At /3 = — 3 we find that out requirement is that the pre- 
correction ellipticity power spectrum must be roughly 1500 
times smaller that the lensing power spectrum at £0 = 1000 
{G'el"" ~ 3 X 10"^° and C^l ~ 1.9 x 10"". This can also be 
expressed as £lCll"'' ~ 3 x 10"" and £lC'/^ ~ 1.9 x 10"''). 

It is, therefore, interesting to consider an extension be- 
yond the simple power-law. We investigate the limits that 
we set on a power spectrum with the following form. 



/32-/3i 



CI 



£0 



+ 1 



to 



(29) 



where /3i is the low £ slope and 132 is the high £ slope. The 
results for this are shown in Figure [H We also note that the 
diagonal that runs from [0,0] to [-4,-4], i.e. Pi = P2 — (3 
corresponds to the power law case shown in Figure [1] 

Not surprisingly, we see that the most relaxed con- 
straints on Gf come from large /3i (less negative) and small 
P2 (more negative). 



5 INTERPRETATION AND CONCLUSIONS 

We have investigated the way that the intrinsic variation 
of the PSF of an observation can be corrected using stars, 
which give a noisy estimate of the PSF at random points 
in the field. This allows us to asses the potential systematic 
contamination that comes from an observation with given 
ellipticity and size powerspectra for the PSF. We show how 
software can be used on large scales to correct the PSF using 
the information form the stars. We have used Wiener filter- 
ing to model the PSF variation. For a Gaussian Random 
field this is known to be the optimal filter for reducing the 
residuals. On small scales however, our ability to calibrate 
the PSF fails due to the finite number of stars available. 
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Figure 1. The upper panel shows a possible example of a PSF 
power spectrum. In this example, we focus on the impact of 
ellipticity variation (explicitly the variation in one of the com- 
ponents of ellipticity). For this example, we assume that both 
components of shear (71 and 72) contribute equally to errors (as 
shown in Equation I26|l . We have also chosen an example where 
the ellipticity variation and size variation contribute equally 
(Fjcr^^^^j = G'^'^syaei)- With these assumptions, the example 
shown here leads to cr'^yc, = 10~^. The red dashed line shows 
the initial ellipticity power spectrum before Wiener filtering with 
/3 = —3. The dotted line shows the noise contributions (C"(£) of 
Eauation ll6l l. The assumptions that go into calculating this noise 
level are given in Section 14.11 The solid red curve shows resid- 
ual power spectrum after filtering. We see that on large scales 
(small £), filtering (i.e. PSF fitting) is able to correct the original 
PSF to a large extent (down to the noise). On small scales (large 
i), the noise becomes dominant. Hence, the filtering process is 
not able to correct PSF, making us dependent on the underly- 
ing, inherent PSF variation. The transition between these two 
regimes - in this example at ^ ~ 500 - marks the point where 
on larger scales the PSF correction depends on the quality of 
the analysis software and on smaller scales PSF variations need 
to be controlled in the observations themselves. The solid black 
curve is an example of the lensing power spectrum (see pub- 
lic code for details). We see that for /3 = —3 our requirement 
is that the pre-correction ellipticity power spectrum be roughly 
1500 times smaller that the lensing power spectrum at £0 = 1000 
C^eT" ~ 3 X 10-1° and ~ 1.9 X 10"", which can also be 
expressed as ^^c^ens ^ 3 ^ ^^^j ^2^ei ^ g ^ iQ-Ty -pj^g 

lower panel shows the post filtering residual power spectrum of 
the errors on the shear using a log-linear linear scale. The integral 



Figure 2. The tolerance on the powerspectra of ellipticity (solid 
black) and radius (red dashed) as a function of the power-law 
slope. We make the same assumptions about the relative strengths 
of each component as we do for Figure 1 . To zeroth order the ratio 
of the powerspectra is given by the ratio F^/G^. We see that 
for both very steep and shallow values of beta, the requirements 
become more stringent. The most relaxed requirements are for 
/3 ~ — 3 (shown in Figure [TJ. 
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component. 



In this regime, lensing measurements become dependent on 
the underlying PSF variation of the instrument. The bur- 
den on these scales then falls on the hardware, since we 
need to perform measurements that have little small scale 
PSF power. This, therefore, provides us with a mechanism 
to make robust predictions about the systematic floor of 
lensing surveys, both current and future. 

To illustrate the power of this approach, we place re- 



Figure 3. The tolerance on ellipticity power spectrum amplitude 
as a function of /3i and [32 , the inner and outer slopes of Equation 
1291 We see that the best results are given for spectra that are 
shallow for low £ and steep for high £. A diagonal cut through this 
figure (bottom left to top right corner) leads to results shown for 
the solid curve of Figure |2l 
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quirements on the size and ellipticity powerspectra that a 
Euchd-like survey would need to stay within the system- 
atic requirem ent of a^ys < 10~^ (jAmara and Refregiedl2008l : 
iKitching et al. 2008). We do this for the case where the er- 
rors on the two components of shear are comparable and 
the contributions to the final systematics from size and el- 
lipticity are comparable. Ideally, the PSF of the observations 
would have a shallow slope on large scales ( /3i > — 3), which 
helps reduce the contributions to sample variance, and a 
steep small scale slope (/32 < —3), which helps contain shot 
noise contributions. For the case where the initial PSF is a 
power- law with f3 = —3 (see Equation I28p . we find that the 
power spectrum of PSF ellipticity at f = 1000 needs to be 
C"^i < 2 X 10~^^ and the power spectrum of size (R^) needs 
tobeC^ < 2xlO"^^This can be compared to the lensing 
power spectrum at this scale of (7'="» ~ 3 x 10"^". Stated 
another way, this means that the pre-correction ellipticity 
power-spectrum needs to be 1500 times smaller than the 
lensing signal (at £ = 1000) for a Euclid-like survey to be 
statics limited rather than systematics limited. In the case 
of space-based surveys, this can be seen as a requirements 
on the auto-correlation power spectrum of the time variable 
instrument PSF. 

For ground-based surveys, if the auto-correlation power 
spectrum of images is found to exceed our requirement, mit- 
igation strategies need to be found since we have only lim- 
ited control over the PSF pattern, which is mostly domi- 
nated by atmospheric effects. For a Gaussian random field 
Winer filtering will give the smallest residuals when mod- 
eling. This makes our predictions optimistic. However for 
a non-Gaussian field other technics , e.g. Principal Compo - 
nent Analysis (PGA) proposed bv Ijarvis and Jalnl l|2008l) , 
may be able to do better than Wiener filtering and reach 
comparable results to Wiener filtering with a Gaussian field. 
This alone may not be sufficient to meet our systematic re- 
quirements since this does not tackle the problem of a high 
initial PSF power spectrum. Other techniques would then 
be needed. As an example of possible options, the shape of 
a galaxy from one image can be correlated with the shape of 
another galaxy in a different image. To set the requirements 
on this we would also need to consider our sensitivity to the 
cross-correlation power spectrum of the PSF between the 
two images. 
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